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A Gaussian fluctuation formula is proved for linear statistics of complex random 
matrices in the case that the statistic is rotationally invariant. For a general linear 
statistic without this symmetry, Coulomb gas theory is used to predict that the 
distribution will again be a Gaussian, with a specific mean and variance. The 
variance splits naturally into a bulk and surface contibution, the latter resulting 
from the long range correlations at the boundary of the support of the eigenvalue 
density. 

The phenomenum of universal conductance fluctuations in mesoscopic wires (see e.g. |]l|]) has 
provided the motivation for a number of theoretical studies into fluctuation formulas for linear 
statistics in random matrix ensembles @. To understand the reason for this, we first recall 
that the striking feature of the conductance fluctuations is that they remain of order unity even 
though the conductance itself is proportional to the number of channels A^. Now, in random 
matrix models of this effect, the conductance can be written as a linear statistic of a certain 
random matrix ensemble (we recall that A is said to be linear statistic of the eigenvalues Xj if 
it can be written in the form A = ^^=ia{Xj) for some random function a). In this setting, 
the theoretical explanation for the phenomenum of universal conductance fluctuations is as an 
example of a universal fluctuation formula in random matrix theory, the first example of which 
was given in the pioneering work of Dyson and Mehta Q . 

For random matrix ensembles in which the support of the density is one-dimensional, for 
example Hermitian or unitary random matrices, (Gaussian) fluctuation formulas are now well 
understood both at a heuristic (see references cited above) and rigorous level |, |, |]. It is the 
purpose of this Letter to initiate the study of fluctuation formulas in complex random matri- 
ces [0, 1^, for which the eigenvalues uniformly fill a disk or ellipse in the complex plane. We 
remark that complex random matrices have occured in recent physical studies of the localization- 
delocalization transition in non-Hermitian quantum mechanics llO] and chiral symmetry break- 



ing in lattice QCD The distribution of a linear statistic is then an observable quantity after 



averaging over many random copies of these systems. 

To begin, we recall [Q] that for a random matrix with complex elements Ujk + ivjk indepen- 
dently distributed with Gaussian distribution ;|e~'^*^l"j'=l^^l^^'=l^\ the corresponding probability 
distribution of the eigenvalues Xj = Xj + iyj is proportional to 
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where fj = {xj,yj). Furthermore, to leading order, the support of the density of the eigenvalues 
is the disk of radius \jN/c. For the purpose of studying fluctuation formulas it is convenient to 
choose c = so that the support of the density is the unit disk. The Fourier transform of the 
distribution of Pr(A = u) is then given by 

nili /r. dn e-m? ni<,<.<Ar 14 - r,f ■ 

Suppose now that the linear statistic is of the form A = Y^^=i '^(l^j Di so that the random 
function a only depends on the distance from the origin. Introducing polar coordinates and 
using the Vandermonde determinant expansion of Y{i<j<k<N{^k — Zj) the integrals in (||) can 
be evaluated with the result 
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This is the exact expression for finite N. To obtain its form for ^ cxd, we change variables 
s — > /s and expand the integrand about its large-Z maximum at s = 1. A straightforward 
calculation then gives 
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with ^ ^ 

^i = 2N [ ra{r)dr, a'^ = - f r{a{r)fdr. (5) 

Thus the distribution of A is a Gaussian with mean and variance as given by (|5|). Note in 
particular that the variance is 0(1). 

Next we address the more general situation in which a is not rotationally invariant. To make 
progress we must proceed heuristically. The p.d.f. (|^) can be interpreted as the Boltzmann factor 
of the two-dimensional one-component plasma (2dOCP) at the special value of the coupling 
r = 2 [^]. Using linear response theory and macroscopic electrostatics, it is possible to argue 
[13, 14] that in general the distribution of a linear statistic in a classical Coulomb system 
in the conductive phase will be Gaussian (this assumes also that the random function varies 
over macroscopic distances relative to the interparticle spacing). The Gaussian distribution is 
uniquely characterized by its mean and variance. But independent of the underlying distribution, 
these quantities are given by 

fi=—[ dfa{r), cr^ = [ dfia{fi) [ dr2a{r2)S{ri,r2), (6) 
vr J A J A J A 

where S{fi,f2) is given in terms of the truncated two particle distribution function by S{fi,f2) = 
pj2^{fi, f2)+ ^^{r—r') , ^ is the particle density and A denotes the unit disk. We see immediately 
from the formula for fi in (^) that the formula for /x in (^) is reclaimed if a(r) = a(|r|). More 
challenging is to reproduce the formula for a"^, and to proceed to generalize this formula for 
general a(r). 
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This task can be undertaken by again appeahng to Coulomb gas theory. In the infinite 
density limit the function S{fi,f2) in for the 2dOCP with general coupling T is expected to 
have the bulk universal form Ifil 
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At r = 2 this can be checked from the p ^ oo limit of the exact formula Q S]^^\]^{fi,f2) = 
— p2e~'^''l''i~'"2p -\- pSi^fi — r2). Substituting (^) in (^), and integrating by parts (ignoring possible 
boundary terms, which are separately treated below) gives 
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In the special case a(r) = a(|f|), F = 2, (^) reproduces the result (|^ for cr^. 

The crucial difference between the case of general a(r) and the rotationally invariant case 
a(r) = is that in the latter case cj^ contains a contribution from the surface correlations 

of the same order (0(1)) as the contribution from the bulk correlations. This effect, due to the 
long-range nature of the correlations at the boundary of Coulomb systems |15], was first noted 
by Choquard |16| and collaborators, who studied the variance of the dipole moment (a(r) = x) 
for classical Coulomb systems. Indeed the variance of this statistic was used to compute from 
microscopic statistical mechanics the macroscopic shape dependent dielectric susceptibility of the 
Coulomb system. 

Like in the bulk, the correlation S{fi,f2) has a universal form for fi and r2 at the surface. 
However, unlike the situation in the bulk, this correlation is long-ranged and shape dependent. 



For A a unit disk, the universal form is |14] 
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where polar coordinates have been introduced. At F 
from the exact evaluation of S'(fi,r2) in the finite system |16]. Substituting in (^) gives 



(9) 



2 this form can be derived explicitly 
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This quantity vanishes for a(r) = a(|r|). 

Consider now a more general ensemble of complex random matrices , in which the mem- 
bers, J say, are of the form J = H + ivA. Here H and A are Gaussian Hermitian random 
matrices with joint p.d.f.'s for the elements proportional to exp ( - ^TiX^^ {X = H,A and 
r = (1 — v'^)/{l + f^)). The corresponding eigenvalue p.d.f. is proportional to 



exp 



^ x2 7/2 , 

j=l ^ l<j<k<N 



(11) 
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(note that in the case r = (11) agrees with (||)), and to leading order the support of the 
eigenvalue density consists of an ellipse with semi-axes A = (1 + r), B = (1 — r). The eigenvalue 
density itself is uniform and thus has the value N/{-k{1 — r^) inside the ellipse. 

The p.d.f. (|Tl| ) can be interpreted as the Boltzmann factor of the 2dOCP at T = 2 in a 



quadrupolar field |17, 18]. Thus, Coulomb gas theory gives that as ^ oo (infinite density 
limit) the distribution of a linear statistic will again be Gaussian. The mean will be given as in 
but with A now the ellipse specifying the support of the eigenvalues, and the factor N/i: 
which represents the eigenvalue density replaced by N/{'ir{l — r^)). With A the ellipse, the bulk 
contribution to the variance is again given by (Mj. For the surface contribution, we require the 



fact that the universal form of the surface correlation at the boundary of an ellipse is [16 

Ssnr.m,Vl),{C2,m) = ^^TTTT^^, slj^^^Og sin ^^^^ ) 6 {^i - 6 {^2 - Cb) , (12) 

T^^^h{(.b,m)h{(b,m)\dmdv2 2 j 

where (^,r/) are elliptic coordinates, specified hy x + iy = cosh(^ + irj), and h{^b,r])dr] gives 
the differential surface element. As in the disk case (^), this form has been explicitly verified 
18| from the known exact expression for S{fi,f2) in the finite system. Since the semi-axes are 



specified by A = cosh^j,, B = sinh.^^, is related to r by tanh^;, = (1 — r)/(l + r). Substituting 
(El) in (I) gives 
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(note the similarity between (13) and (|To|)). 

There is a simple linear statistic for which the exact distribution can be calculated, thus 
allowing the above predictions to be tested. This statistic is the linear function a{x, y) = 
ciqx + coiy. Substituting in the analogue of (^) for the p.d.f. (pT|) raised to the power r/2, 
the resulting dependence on k can be simply calculated by completing the square and changing 
variables in the integrand (see for an analogous calculation in the case of Hermitian random 
matrices). We find 

independent of N. Thus cr^ = ^(cfo(l + ''") + Coi(l — t)). Substituting the linear function a{x,y) 



in (g) and ( p!3|) and adding the result verifies that the general formulas reproduce the exact 
value. 
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